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ABSTRACT 


Comparisons  between  the  CO^  telluric  absorptions  observed 
in  high  altitude  solar  spectra  and  those  predicted  on  a  theoret¬ 
ical  basis  indicates  that  the  atmospheric  concentration  of 
decreases  above  the  tropopause.  CO^  mixing  ratios  were  deter¬ 
mined  on  the  basis  of  the  absorptions  observed  in  the  2.7u  and 

4.3m,  spectral  regions.  Ihe  concentration  determined  by  oth  bands 

-4 

gives  a  value  of  2.1-10  as  the  volume  concentration  of  above 

-4 

17  km  as  compared  with  the  value  of  3.3*10  for  the  surface  con¬ 


centration  . 
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1,  INTRODUCTION 

It  has  long  been  assumed  that  the  concentration  of  carbon 
dioxide  in  the  atmosphere  is  independent  of  altitude  at  least 
through  the  stratosphere  and  varies  only  slightly  from  its  mean 
value  with  location  and  time.  A  detailed  analysis  of  infrared 
solar  spectra  obtained  at  various  altitudes  by  means  of  a  balloon 
borne  spectrometer  system  indicates  that  there  is  a  large  decrease 
in  concentration  of  carbon  dioxide  above  the  tropopause. 

1  2 

In  previous  studies  *  the  atmospheric  absorptions  due  to 
CO2  in  the  2.7ji  anci  4.3|j,  regions  observed  at  various  altitudes  were 
compared  with  the  theoretical  absorptions  calculated  using  a  number 
of  techniques  without  arriving  at  any  conclusions  concerning  the 
concentration  of  002*  Since  these  earlier  publications  a  number 
of  advances  have  been  made  in  computers,  computational  procedures 
and  laboratory  data  concerning  molecular  absorption.  Additional 
experimental  data  have  also  been  obtained.  Using  the  new  techniques 
and  the  new  data  it  has  been  possible  to  determine  the  concentration 
of  CO2  in  the  lower  stratosphere  from  the  observed  atmospheric  absorp¬ 
tions  with  sufficient  accuracy  to  show  that  the  concentration  decreases 
significantly  above  the  tropopause. 

2.  EXPERIMENTAL  DATA 

The  high  altitude  solar  spectra  used  in  this  study  were 
obtained  by  means  of  a  balloon  borne  spectrometer  system  which  is 
flown  in  conjunction  with  a  biaxial  pointing  control  which  maintains 


2 


the  solar  radiation  on  the  apectrometer  entrance  slit.  The  instru- 

3 

mentation  has  been  described  in  detail  in  a  previous  publication. 

The  absorption  bands  used  in  this  study  were  the  two  vibration  ro- 
tation  bands  in  the  vicinity  of  2.7^  and  the  strong  absorption  at 
A.3u.  The  data  concerning  the  variation  of  the  2.7u  absorption  in 
the  solar  spectrum  with  altitude  above  the  earth's  surface  were  ob¬ 
tained  during  a  balloon  flight  made  September  30,  1965  from  Chico 
California.  Data  concerning  the  variation  of  the  4.3^  band  with 
altitude  were  obtained  during  a  balloon  flight  made  May  14,  1964 
from  Holloman  AFB,  New  Mexico.  During  both  flights  solar  spectra 
were  recorded  at  various  altitudes  from  the  ground  through  30  km. 

3.  ANALYSIS  OF  THE  ABSORPTIONS 

The  2 . 7p.  band  of  CD^  is  overlapped  by  water  vapor  absorption 
bands  at  the  same  wavelength,  but  the  decreased  amount  of  water 
vapor  in  the  atmosphere  above  the  tropopause  allows  the  overlapping 
to  be  taken  into  account.  At  levels  below  11  km  however,  the  water 
vapor  absorptions  are  large  enough  to  make  it  impractical  to  deter¬ 
mine  the  carbon  dioxide  concentration  using  this  absorption  region. 

At  the  higher  altitudes  the  4.3p,  band  is  not  affected  by  the  absorp¬ 
tion  of  other  atmospheric  gases.  Only  the  region  from  2300  -  2400  cm 
was  used  in  the  determination  of  the  00£  concentration  since  it  is 
felt  that  the  smaller  contributions  of  isotopic  and  "hot"  states 
in  this  region  would  result  in  a  more  accurate  determination  of  the 
00^  concentration.  The  great  intensity  of  this  band  causes  the 
absorption  to  be  almost  total,  at  least  over  a  large  part  of  the 
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wavelength  interval,  for  altitudes  less  than  23  km.  In  view  of 
these  limitations  the  mixing  ratio  was  determined  only  for  alti¬ 
tudes  above  11  km,  and  by  use  of  two  different  bands  above  23  km. 

The  infrared  absorptions  considered  here  depend  on  the 
pressure  environment  and  temperature  of  the  absorbing  gas.  In 
order  to  predict  theoretically  the  absorptions  to  be  expected 
under  atmospheric  conditions  it  is  necessary  to  incorporate  some 
techniques  for  treating  the  "slant  path"  calculations.  The  var¬ 
iation  in  pressure  environment  of  the  gas  was  treated  by  breaking 
up  the  atmosphere  into  incremental  layers.  In  each  layer  the 
pressure  was  assumed  to  be  constant  and  the  transmission  was  cal¬ 
culated  for  each  layer.  The  temperature  dependence  of  the  absorp¬ 
tion  arises  from  the  dependence  of  the  line  strength  and  half-width 
on  temperature.  The  overall  dependence  of  the  integrated  absorb- 
tance  on  temperature  is  at  most  as  T^  for  the  bands  considered 
here.  In  view  of  the  slow  dependence  on  temperature  the  atmos¬ 
pheric  temperature  was  taken  to  be  constant  at  236°K.  This 
temperature  is  within  25°K  of  the  atmospheric  temperature  for 
at  least  957.  of  the  absorber  concentration.  In  the  upper  strato¬ 
sphere  the  temperature  reaches  higher  values,  however,  the  contri¬ 
bution  of  the  00 2  in  this  region  to  the  absorption  is  small.  The 
error  introduced  by  this  constant  temperature  assumption  is  small 
and  is  a  maximum  at  the  lower  altitudes.  At  30  km  the  error 
introduced  should  be  almost  zero. 

The  absorption  parameters  -  line  frequencies,  strengths, 
and  energies  for  the  4.3  micron  band  were  generated  from  the 
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rotational  constants,  band  centers  and  band  strengths  published  by 

A 

Gray.  Ihe  only  modification  of  these  values  was  the  shifting  of 

1 2  18 

the  band  center  of  the  fundamental  isotopic  state  C  Oj  from 
2333.00  cm  ^  to  2332.16  cm  *  on  the  basis  of  a  recent  paper. ^  Ihe 
same  parameters  for  the  2.7m.  region  were  generated  from  the  band 
constants  given  by  Calfee  and  Benedict.^  The  values  were  generated 
rather  than  using  the  listed  values.  Some  errors  in  line  strengths 
presented  by  these  authors  seem  to  be  due  to  a  computer  malfunction. 

This  is  evident  in  the  fact  that  in  some  bands  the  sum  of  strengths 
of  the  lines  is  not  equal  to  the  listed  band  strengths. 

The  calculations  were  performed  on  a  line  by  line  basis. 

All  the  lines  in  the  bands  were  taken  into  account,  and  an  extinction 
coefficient  determined  at  points  0.001  cm  apart.  This  extinction  co¬ 
efficient  took  into  account  the  Doppler  effects  on  the  Lorentz  lines. 
Once  the  extinction  coefficient  was  known,  the  absorption  at  each 
point  was  determined.  After  degrading  this  monochromatic  absorption 
by  the  application  of  a  triangular  slit  function,  the  integrated  ab¬ 
sorption  was  determined.  The  application  of  the  slit  function  reduced 
the  errors  which  could  have  occurred  at  the  ends  of  the  interval  of 
integration  were  it  not  used. 

The  accuracy  of  the  method  of  calculating  the  absorption  to 
be  expected  and  various  absorption  parameters  was  checked  by  comparing 
the  laboratory  spectra  given  by  Burch^  with  those  calculated  for  similar 
paths  and  pressures.  The  result  of  the  comparison  is  given  in  Table  I. 
Better  agreement  is  obtained  by  taking  the  half- widths  of  the  lines  to 
be  .07  cm  ^  at  one  atmosphere  pressure  Instead  of  .08  cm  ^  as  suggested 
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by  Calfee  and  Benedict.  Because  of  the  superior  agreement  achieved 
in  the  3650-3766  cm  ^  interval,  a  subinterval  in  this  region,  3664- 
3766  cm  *  was  used  for  the  determination  of  the  atmospheric  CO2  con¬ 
centration.  It  should  be  noted  also  that  the  use  of  a  too  small 
value  for  the  half- width  of  the  line  would  result  in  too  large  rather 
than  too  small  a  mixing  ratio  for  the  CO2. 

The  water  vapor  absorption  in  the  2.7m*  region  was  taken  into 
account  by  what  might  be  described  as  a  reversal  of  the  Curtis- Godson 
approximation.  In  normal  use  of  the  Curtis- Godson  approximation,  the 
amount  of  absorbing  material  is  known  and  a  scaled  pressure  is  used 
as  though  the  absorption  was  due  to  a  constant  pressure  path  through 
the  given  amount  of  absorber.  In  the  present  work  the  integrated 
absorption  in  the  interval  3764  cm  *  to  3782  cm  \  which  is  entirely 
due  to  H^O,  was  used  to  determine  the  amount  of  water  vapor  which 
would  have  to  be  present  at  30  km  altitude  to  account  for  the  observed 
value  of  the  integrated  absorption.  Thus  rather  than  scaling  the 
pressure  with  a  constant  amount  of  absorber  the  pressure  was  kept 
constant  and  the  quantity  of  absorber  was  scaled.  No  physical  sig¬ 
nificance  can  be  attributed  to  the  scaled  quantity  of  absorber,  but 
it  can  be  used  to  predict  the  quantity  of  water  vapor  absorption  in 
the  region  of  absorption  by  carbon  dioxide.  Its  validity  rests  on 
the  same  arguments  as  the  Cur tis-Godson  approximation.  It  is  im¬ 
portant  to  keep  in  mind  that  this  procedure  is  only  being  used  at 
altitudes  where  the  absorption  by  water  vapor  alone  is  much  less  than 
that  due  to  carbon  dioxide  alone.  No  such  procedure  is  required  for 
the  4.3p,  band  since  it  is  not  overlapped  by  other  bands  ir  the  region 


of  interest  here. 


The  slant  path  calculation  was  performed  as  follows:  The 
atmosphere  above  a  given  altitude  at  which  the  pressure  is  p  is 
made  of  ten  constant  pressure  layers  at  pressures  O.lp,  0.2p  •••• 
••••0.9p,  and  l.Op.  The  mixing  ratio  was  taken  as  the  same  in  all 
layers.  The  lines  in  the  2.7y,  region  were  taken  to  have  a  pure 
Lorentz  shape  in  each  layer  and  this  shape  was  then  convoluted  with 
the  Gaussian  shape  of  a  pure  Doppler  broadened  line.  The  product 
of  the  resulting  line  shape  and  the  absorber  concentration  was  summed 
over  the  ten  levels.  The  same  procedure  was  used  for  the  4.3|jl  band 
except  a  modified  rather  than  pure  Lorentz  line  shape  was  used.  The 


stant3  used  were  a*0.46  and  b®0,46,  A  line  shape  of  this  type 
would  be  more  accurate  for  the  2.7m.  band  but  the  values  of  the 
constants  a  and  b  are  not  known  for  this  band.  Calculations  have 
indicated  that  the  error  caused  by  this  omission  of  the  modifica¬ 
tion  in  the  line  shape  will  be  small.  The  errors  will  be  greatest 
at  low  altitudes  where  the  contribution  of  the  wings  of  the  lines 
is  more  important.  The  agreement  obtained  between  the  two  bands 
at  the  higher  altitudes  confirms  the  validity  of  use  of  the  pure 
Lorentz  line  shape  convoluted  with  the  Doppler  shape  for  the 
upper  altitudes.  It  should  again  be  emphasized  that  all  errors 
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in  line  shape  (including  the  possibility  of  an  incorrect  value  for 
the  half-width)  result  in  too  large  a  mixing  ratio  rather  than  too 
small . 

The  primary  results  of  the  calculations  are  given  in  Table 
II.  The  integrated  absorption  between  3664  cm  *  and  3764  cm  *  for 
the  2.7}i  band  and  between  2300  cm  *  and  2400  cm  ^  for  the  4.3p»  band 
is  due  to  the  absorptance  along  a  path  making  an  angle  0  with  zenith 
from  the  given  altitude.  The  mixing  ratio  was  determined  by  perform¬ 
ing  the  calculations  described  above  for  different  values  of  the 
mixing  ratio  and  then  interpolating  between  the  values  to  determine 
the  mixing  ratio  which  would  give  agreement  with  the  observed  value 
of  the  integrated  absorptance.  The  variation  of  mixing  ratio  with 
altitude  can  be  seen  better  in  Figure  1  which  is  a  plot  of  effective 
mixing  ratio  vs  altitude.  The  lowest  altitude  for  which  an  effective 
mixing  ratio  was  determined  is  11.87  km  where  the  mixing  ratio  has 
already  decreased  by  10%.  However  it  should  be  emphasized  that  the 
effective  mixing  ratio  is  the  mixing  ratio  which  if  applied  through¬ 
out  the  atmosphere  will  give  the  observed  value  of  the  integrated 

-4 

absorption.  Thus  at  11.87  km  an  effective  mixing  ratio  of  2.74*10 

fractional  concentration  by  volume  is  found  but  at  altitudes  above 

-4 

17  km  a  value  of  2.15*10  is  found.  This  means  that  the  actual 

-4 

mixing  ratio  at  11.87  km  is  really  greater  than  2.74*10  since  too 
large  a  value  of  the  mixing  ratio  has  been  used  in  the  upper  layers 
of  the  calculation  for  the  11.87  km  case. 

If  the  simplifying  assumption  is  made  that  the  absorption  A 
is  a  function  of  the  integral  of  the  atmospheric  pressure  p  with 
respect  to  the  optical  path  u,  that  is 

A  =  A(R)  where  R  =  pdu 
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then  the  mixing  ratio  versus  altitude  profile  can  be  approximated. 
The  mixing  ratio  obtained  by  assuming  that  the  COj  is  uniformly 
distributed  above  17  km  and  modifying  the  mixing  ratio  at  lower 
altitudes  so  as  to  give  the  same  value  of  R  as  used  in  the  spectral 
calculations  is  shown  as  the  dashed  line  in  Figure  1.  It  can  be 
seen  that  the  variation  in  mixing  ratio  occurs  very  close  to  the 
tropopause. 

4.  DISCUSSION 

Most  calculations  of  slant  path  transmissions  in  the  past 
have  been  done  on  the  basis  of  the  Curtis-Godson  approximation. 

In  order  to  see  how  this  approximation  would  affect  the  results 
obtained,  the  integrated  absorptance  in  the  2.7p,  region  was  cal¬ 
culated  both  by  the  use  of  the  Curtis-Godson  approximation  and  by 

the  use  of  the  ten  layer  model  with  Doppler  broadening  in  each 

-4 

level  and  a  mixing  ratio  of  3.0*10  in  all  cases.  The  integra¬ 
ted  absorptances  obtained  with  the  Curtifr-Godson  approximation 
were  too  large  at  the  lower  altitudes  and  too  small  at  the  upper 
altitudes.  The  results  are  reviewed  in  Table  III.  The  change 
required  in  the  effective  mixing  ratio  to  bring  the  Curtis-Godson 
approximation  into  agreement  with  the  experimental  data  is  also 
shown*  The  differences  at  the  lower  altitudes  are  independent 
of  the  Doppler  effect  and  are  due  to  insufficiencies  in  the 
approximation  while  those  at  the  upper  altitudes  are  mainly  due 
to  Doppler  effects. 

No  explanation  for  the  decrease  in  the  mixing  ratio  of 
carbon  dioxide  can  be  offered  at  this  time.  Although  the 
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molecular  weight  of  CO^  is  large  the  available  information  concern¬ 
ing  gravitational  separation  indicates  that  the  altitude  at  which 
the  decrease  in  mixing  ratio  occurs  is  much  too  low  to  be  due  to 
this  effect.  Tile  occurrence  of  the  decrease  at  the  tropopause 
seems  to  be  significant  and  appears  to  indicate  that  the  amount  of 
COj  transported  through  the  tropopause  is  small,  or  there  is  some 
reaction  which  destroys  carbon  dioxide  in  the  lower  stratosphere. 


Figure  1.  Effective  CX) 2  Mixing  ratio  as  determined  from  atmospheric  absorptions  as  observed  at  various 
altitudes.  Procedure  for  determining  the  refined  mixing  ratio  is  discussed  in  the  text. 


— _ 
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Table  I 


Interval 

Pa  th 
Atmos 

CM 

Equivalent 
Pressure 
mm  Hg 

Integrated  Absorptance  cm  ^ 
Experimental  Calculated 

Half-width  at 
Atmos  cm 

i 

50.33 

77.52 

67.02 

68.05 

0.07 

i 

50.33 

77.52 

67.02 

70.93 

0.08 

i 

50.29 

762.05 

104.23 

108.00 

0.08 

i 

16.91 

762.05 

78.77 

84.62 

0.08 

i 

0. 3047 

817.30 

7.16 

9.07 

0.07 

i 

0.3047 

817.30 

7.16 

9.34 

0.08 

2 

50.33 

77.52 

63.70 

61.73 

0.07 

2 

50.33 

77.52 

63.70 

64.03 

0.08 

2 

50.29 

762.05 

88.42 

91.17 

0.08 

2 

16.91 

762.05 

74.22 

75.43 

0.08 

2 

0.3047 

817.30 

10.38 

10.45 

0.07 

2 

0.3047 

817.30 

10.38 

10.83 

0.08 

Interval  1: 

3490.3  - 

3650.8 

-  1 
cm 

Interval  2: 

3650.8  - 

3766.0 

-1 

cm 
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Table  II 


Band 

Altitude 

km 

Sec  6 

.  -1 
Avcm 

Mixing  Ratio  x  10 

2.7 

11.9 

1.43 

70.4 

2.74 

2.7 

13.6 

1.45 

61.5 

2.56 

2.7 

16.9 

1.49 

42.5 

2.14 

2.7 

18.0 

1.51 

35.5 

2.29 

2.7 

20.3 

1.56 

27.0 

2.04 

2.7 

22.2 

1.59 

23.2 

1.98 

4.3 

26.6 

1.78 

56.0 

2.22 

4.3 

27.4 

1.73 

52.0 

2.16 

4.3 

27.9 

1.70 

50.0 

2.16 

4.3 

29.2 

1.63 

43.0 

2.01 

2.7 

29.7 

1.88 

9.5 

2.16 

4.3 

30.4 

1.53 

37.0 

1.83 
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Table 


A  cm 
v 

Altitude(km)  Curtis- Godson 


11.87 

70.68 

13.56 

63.00 

16.87 

45. 36 

18.00 

39.56 

20.32 

30.59 

22.23 

25.02 

III 


Av  cm 

10  Level 

Difference  in 

4 

Mixing  Ratio  x  10 

70.21 

-0.10 

62.53 

-0.09 

45.11 

-0.04 

39.40 

-0.03 

30.69 

0.01 

25.24 

0.05 
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